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Abstract 

We investigate the multi-soliton solutions to tlie generalized discrete KdV equa- 
tion. In some cases a soliton with smaller amplitude moves faster than that with 
larger amplitude unlike the soliton solutions of the KdV equation. This phe- 
nomenon is intuitively understood from its ultradiscrete limit, where the system 
turns to the box ball system with a carrier. 
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1 Introduction 

The discrete KdV equation 
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is an integrable partial difference equationJT]. Here n, f e Z and (5 e R is a parameter. 
The eq. ([T]i turns to the continuous KdV equation by taking an appropriate continuous 
limit and has multi-soliton solutions. The solitons in ([T]i have similar properties to 
those in continuous KdV equation. They do not change their amplitude after collision 
and a soliton with larger amplitude moves faster than that with smaller amplitude. Note 
that when we put 

we obtain equivalent coupled equations 
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In 1990, Takahashi and Satsuma fl] proposed a cellular automaton with soliton 
solutions which is now called a 'box ball system' (BBS). The time evolution rule of 
BBS is as follows: 
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where U'„ is a number of balls in «th box at time f and Cb means every box holds at 
most Cb balls. The parameter Cb is called 'capacity of the box'. We introduce another 
variable to obtain coupled equations which are equivalent to the BBS O): 



This system is known to be derived from the discrete KdV eq. ([U by a limiting pro- 
cedure called the ultradiscretization|3|. The BBS also has multi-soliton solutions. A 
peculiar feature of solitons in BBS is that the velocity of a soliton is exactly propor- 
tional to its amplitude, the number of balls constituting the soliton, just like a soliton in 
the continuous KdV equations. 

The 'box ball system with a carrier' (BBSC) [4| is a generalization of the BBS 
defined as follows: 



where Cc is a positive parameter which is called 'capacity of the carrier'. To describe 
the evolution rule of the BBSC we prepare a 'carrier' of balls. We assume that the 
carrier can carry at most Cc balls. From time f to f + 1, carrier moves from the -oo site 
to the oo site and passes each box from the left to the right. The carrier gets as many 
balls as possible from the nth box and, at the same time puts as many balls as possible 
to the vacant space of the n + 1th box for n = ■ ■ ■ , -2, -1, 0, 1, 2, ■ ■ ■ . The action of the 
carrier and the time evolution of BBSC for Cb - 3 and Cc = 4 is described in fig. [T] 
and fig. |2] 

In this paper we treat a generalized system of (|2]) and calculate the velocity and the 
amplitude of the soliton solution to the equation. We present several examples of these 
solutions. The ultradiscrete Umit of the generalized equation we treat in this article is 
proved to be equivalent to the BBSC. One interesting property of the solitons in the 
generalized system is that a soliton with smaller amplitude moves faster than that with 
larger amplitude in some parameter region. This phenomenon is intuitively understood 
by the corresponding solitons in the BBSC. 
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min [Cb - U'„, V'„] + max [0, U'„ + V'„- Cc] , 



(5) 



2 



XX 



C„=3 



Cc = 4 



Figure 1 : The action of the carrier exchanging the balls with the box. From time t to 
f + 1, the carrier passes each box from the left to the right. Note that getting the balls 
from the box and putting the balls to the box take place at the same time. 
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Figure 2: The time evolution of the BBSC for Cb = 3 and Cc - 4. 
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2 Generalized discrete KdV equation and its A/^-soliton 
solutions 

Let us consider a generalized discrete KdV equation||5l 



- a) + ax„y„ 
{I - a) + ax'„y'„ , 
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By a scaling transformation 



and change of variables a - a( l - /3) and b 
Yang-Baxter map. Tab ■ (m, v) — > (u', v') 



■Pyl (7) 

- /3(l - a), the eq. ([Sjl becomes the 
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By putting 
we obtain 
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and by taking the limit b — > +oo, 
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Therefore the eq. (|6]l is indeed a generalization of the discrete KdV eq. (|2|. The eq. 
(|6]l is known to be derived from the reduction of the discrete KP equation. 



Proposition 1 (Date-Jimbo-MiwaQ) 

Let us consider 4 component discrete KP equation: 

(fll - b)Tl,tTn +(b- c)TlJ,n + (c - fll)Ti|„T, = 0, 

(ai - b)ThtT„ + {b- c)t;jT,„ + (c - fl2)T;,„T, = 0. 



(9) 
(10) 



Here t - T(li,l2,t,n) {{li,l2,t,n) e Z"*) is the T-function, and ai, 02, b, c are ar- 
bitrary parameters and we use the abbreviated form, r = T(li,l2,t,n), ti^ = t(/i + 
1, /2, t, n), Ti^t = t(Ii + l,l2,t + 1, «) and so on. Then the N-soliton solution to ^ and 
(doll is 



det 

l<iJ<N 



qt - a2 



Pi- qjXPi- <^il \Pi-a2l \pi-bj \pi-c 



(11) 



where {pi,qi]^^i are the parameters distinct from each other and {ji]^^^ are arbitrary 
parameters. 
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Using this fact we obtain the A^-soHton solutions to © by a similar reduction adopted 
in Kakei-Nimmo-Willox |6|. 

Proposition 2 (Kanki-Mada-Tokihirod)) 

The N-soUton solutions of ^ are 



fgn 
gfn ' 



t _ sft_ 
fgt 



where 



det 

l</,i<A' 



- det 

i</,;<A' 



5ij + 



5ij + 



-Pi +15 



Pj + A\pi + I - a / \ -pi + a 



1-/3 



-A- Pi I -Pi+P \ I Pi + 1-/3 



Pi + Pj + A Pi \pi + l-aj \ -pi + a 
with A - I - a - /3. 

Sketch of the proof Imposing the reduction condition 

T/,/, = T, 

to © and ([Tol l gives the constraint 

(ai - pi\la2- Pi , _ ^ 



(12) 

(13) 

(14) 



(15) 



01 - qi/\a2- 

to the parameters [pi, qi] in (fTTl i. Since pi + the constraint becomes /?,■ + qi = 
a\ +02- Then we define / := t, ^ := t;, from the solution to (|9]l and (fTOl i. Putting 

c — a\ 02 — b . [ J— , 

a :- ;-, /3 :— , we find that jcj,, yj, defined by il2\ satisfy the eq. (|6]l. By 



c-b 



c - b 



redefining the parameters as 
the result. 



Pi-ai 
c-b 



c-b 



Ji and by putting li = Z2 we have 

□ 



Let us consider the one-soliton solution (fTSl l and (fl4l) for - I, p\ - p and y\ - y. If 
we take 



A = 



-P+/3 
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, D = 



-A-p 



p + l-a —p + a 2p + A 

the solution of the eq. (|6]l has the following form: 

, _ (1 + CA'B")(1 + CDA'B"^^) 
^" ~ (1 +CDA'B")(1 +CA'B«+i)' 

We assume that the parameters satisfy the condition 

0<Q'<1, 0<y8<l,0</?<Q'+yS-l, 0<r (75- 



a+/3-l' 
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to assure that we have A>0, B>0, C>0 and D > 0, which are sufficient to obtain a 
one-soHton with a bounded ampHtude. The velocity v(p) of x'„ is 



v(p) :-- 



log A 



logB 
1 P 



P + 
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(16) 



We then calculate the fluctuation of with respect to X := B" . With some parallel 
displacement with respect to t we can assume that CA' — 1 without changing the 
amplitude of the solution. We obtain 



dx\. 



1 - BDX^ 



(1 + Bxy(i + Dxy 



Note that 



(B-l)(D-l) 



Therefore x'^ has a local minimum at X = 
A 



<0 p*--, 
if p ^ , and x' 

Ibi5 2 



1 for all X if 



-. Summing up the results obtained above, the amplitude W{p) of one-soliton 



x'„ is as follows: 



W{p) = \x' - lU 



(17) 



where no satisfies B" 



BD 



Theorem 1 

For the velocity v(p) and the amplitude W{p) of one-soliton solution of the generalized 
KdV eq. ^ we have the following properties: 

• W{p) is monotone decreasing on [0, {a + /3 - l)/2] and monotone increasing on 
[(a+/3-l)/2,a+/3- 1]. 

• IfQ<a</3<l then v{p) is monotone increasing on [Q,{a + /3 - l)/2] and 
monotone decreasing on [(a + p - l)/2, a + /3 - 1]. 

• /fO<y6<Qf<l then v{p) is monotone decreasing on [0, (a + p - l)/2] and 
monotone increasing on {{a +/3 - l)/2,a + /3 - I]. 



IfO < a - p < I then v(p) is constantly equal to 1. 
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Sketch of the Proof Basically we just have to consider the fluctuations of v(p) and 
W(p) using elementary calculus. □ 



The following corollary follows immediately. 

Corollary 1 

Wfe have 

. v(p) I v(q) ^ W(p) I W(q) ( if a > P), 

. v{p) I v{q) ^ W{p) I W{q) ( ifa<l3). 

between the velocity and the amplitude of one-soliton solution of (|6|l. 

Therefore if a < /3 the smaller soliton moves faster than the larger one in A^-soliton 
solutions. For example, if we consider the two-soliton solution with parameters a - 
i y3 = B' /'I " ^' ^1 " -5' ^2 = ^, 72 = -35, then we obtain 

log I log I 

v(pi) = = 0.783 > y(p2) = = 0.723, 
log I log 8 

W{pi) = 0.363 < W{p2) = 0.722. 

Fig. |3] shows the two-soliton solution described above. Lastly we note on the case 
a = yS. In this case, clearly every solution has a speed 1 everywhere, since (|6]l becomes 
just xj,^' = y'lj and yj^^j - xj,. In fig. |4]we show a two-soliton solution where the 
parameters are taken as a - f5 - ^, pi - j^, yi - -20, p2 = T2 = 



3 Ultradiscrete limit of the generalized discrete KdV 
equation 

Next we take the ultradiscrete limit of the generalized discrete KdV eq. ^ and the 
A^-soliton solutions. We show that the ultradiscrete limit of ^ gives a box ball system 
with a carrier (BBSC). 

Proposition 3 

The ultradiscretization of Q is equivalent to BBSC ©. 
Proof We define X;,, F,;, A, Bby 




a - exp^-— j, y8 = exp^- — j. 



where e > is an arbitrary parameter, and then take the ultradiscrete limit of (|6]l, which 
means that we take lim elog( ) of both sides of the eq. (|6]l. From the first eq. of (|6]l 
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Figure 3: The two-soliton solution of the generalized discrete KdV equation where 
a - ^ < P - ^. We observe that the smaller soliton takes over the larger one. 
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Figure 4: The solution of the generalized discrete KdV equation where a - /3 = |. 
Every point moves at speed one if a - f5. 
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we obtain 

X'^^ = min [O, B + X'„ + F,',] + Y'„- min [o, A + X,', + Y',^ 
= Y'„ + min [O, B+X'„ + Y'„\ + max [o, -A - X'„ - Y'^ 
= min [-X'„, B + Y',] +X'„ + max [f^, -A - Z^] 
= min[-X,;, B + y;,] + max [x'„ + Y'„ + A, o] - A. 

Note that since < ff < 1, we have A > and Hm^^+o f log (l -exp^-^)) - 
follows. The same fact applies to /3. The relation x'^^y'^^^ - x'^y^ shows that F^^j - 
X'„ + Y'„- . If we displace x;, and Y'^sls " 

U'„^X'„+A, V'„^Y'„ + B, (18) 

we obtain the following equation 

Ui;^ = min [A - U'„, V'„] + max [0, U'„ + V'„ - B] , 
V = U' +V' - U'^^ 

n+1 n ^ ' n ^« ' 

which is exactly the BBSC (EJ itself with Cb- A and Cc - B. □ 

Note that transformations ^ between the discrete systems correspond to ( fTST l in the 
ultradiscrete case and that the ultradiscrete limit of the Yang-Baxter map (|8]l is the 
original BBSC ©. 

Theorem 2 

We have the relation 

l3%a^CB = Cc 

between the parameters of the generalized discrete KdV eq. ^ and those of the BBSC 
^ when < a, /3 < I. 

Proof ff < jS is equivalent to Cc = B < A = Cb- The rest of the theorem is shown in 
the same manner □ 

Thus the inequality relation between the parameters a and /3 of generalized discrete 
KdV equation corresponds exactly to that of the parameters Cb and Cc of BBSC. The 
following example illustrates the soliton solutions of the BBSC when Cc < Cb- The 
fig. |5] shows the time evolution of U'„ if Cb - 3 and Cc - I- The speed of the soliton 

on the left (larger) is - and that of the soliton on the right (smaller) is 1 . The reason 

why the smaller soliton moves faster is clear; the capacity of the carrier is just one and 
it can carry a soliton with amplitude 1 at a time. For a soliton with amplitude 3, the 
carrier cannot carry it all at once and it takes three time steps to move them to the next 
box. 

We conclude that these phenomena in soliton solutions of the BBSC coincide with 
those in the generalized discrete KdV equation described in fig. |3] 
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Figure 5: The time evolution of the BBSC for Cb = 3 and Cc = 1. 
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4 Concluding remarks 



We have obtained the velocity and the amplitude of soliton solutions to the generalised 
discrete KdV equation. We have found the cases in which the smaller solitons take 
over the larger ones. Through the ultradiscrete limit, these soliton solutions turn to 
the solutions of the BBSC. In BBSC, larger cluster of balls can move slower than the 
smaller ones, which corresponds to the phenomena we have obtained for the discrete 
equation in this article. Detailed analysis of the solutions of other discrete integrable 
equations and of its ultradiscrete limits are the problems we would like to address in 
the future. 
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